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PREFACE 


This  volume  includes  papers  presented  during  the  XXI Xth  Congress  and  Plenary 
Assembly  of  the  International  Commission  for  the  Scientific  Study  of  the 
Mediterranean  Sea  (CIESM)  in  Lucerne,  Switzerland,  October  11-19,  1984.  The 
authors  are  a  loose  association  of  scientists  from  the  Donde  V a?  project  and  from 
NORDA.  The  papers  will  also  be  reported  in  the  Rapports  et  Process-Verbaux  of 
CIESM,  where  they  will  be  condensed  to  meet  the  prescribed  two-page  limit  of  that 
report.  Their  format  here  is  enlarged  primarily  by  the  inclusion  of  the  figures 
that  accompanied  the  original  presentations. 

Recent  activities  in  the  Mediterranean  included  the  international  Donde  Va? 
experiment  in  1982,  which  focused  on  the  Alboran  Sea  gyre.  The  extensive  data 
collected  at  that  time  have  now  undergone  preliminary  processing  and  form  the 
basis  for  a  number  of  papers,  including  the  first  five  in  the  present  collection. 
Further  exploitation  of  this  data  is  still  underway,  as  is  development  of  the 
numerical  studies  associated  with  it. 

A  number  of  studies  are  presently  planned  for  the  Mediterranean.  Observational 
efforts  will  reach  a  peak  in  1986,  and  analysis  and  modeling  will  likely  continue 
through  the  end  of  the  decade.  The  total  list  of  planned  work  is  impressive: 
multinational,  multidisciplinary  programs  are  scheduled  in  the  Eastern  and  in  the 
Western  Mediterranean  Basins;  a  comprehensive  study  of  the  Strait  of  Gibraltar 
will  be  supported  by  ONR;  and  continuing  NORDA  programs  (Chemical  Fronts,  Fine 
Scale  Variability)  anticipate  component  programs  in  the  Western  Mediterranean. 

The  final  three  papers  presented  here  are  the  first  of  NORDA' s  contributions  to 
these  programs. 

Present  emphasis  on  the  Mediterranean  reflects  its  importance  in  international 
affairs  and  its  value  as  a  natural  laboratory  where  the  same  physical  processes 
that  control  the  large  ocean  basins  may  be  conveniently  studied.  The  present 
volume  is  thus  a  small  step  toward  a  grand  objective.  In  this  spirit,  we  submit 
it  to  our  colleagues  and  our  sponsors. 


The  Editor 
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The  Donde  V a  Group 

This  article  appeared  in 
the  Ocean  Report,  EOS,  65,  682-684 
September  4,  1984 

INTRODUCTION 

During  June-October  1982  an  international  consortium  of  oceanographers  (see  Table) 
studied  the  circulation  of  the  western  Alboran  Sea  in  an  experiment  entitled 
"Donde  Va?" .  Although  the  English  translation  of  this  title  is  "Where  does  it 
go?",  our  goals  were  more  ambitious  than  the  name  implies.  In  this  overview  of 
Donde  Va?  we  will  discuss  the  oceanographic  background  of  the  region,  the 
objectives  of  the  experiment,  and  the  preliminary  results. 

OCEANOGRAPHIC  BACKGROUND 

The  Alboran  Sea  is  the  westernmost  of  the  many  basins  that  comprise  the 
Mediterranean  Sea.  A  narrow  (10-30  km)  current  of  Atlantic  Water  (fresher  than 
36.5  salinity)  flows  eastward  through  the  Strait  of  Gibraltar  with  a  volume 
transport  of  about  1.4  x  10  m  /sec,  while  saline  Mediterranean  Water  (38.4 
salinity)  flows  westward  beneath  it  with  a  transport  that  is  about  4%  less.  This 
two-layer  flow  maintains  the  salt  and  water  balance  of  the  Mediterranean  which 
annually  loses  about  1  m  over  its  surface  through  the  excess  of  evaporation  over 
river  runoff  and  precipitation  (Lacombe,  1984).  While  intuition  might  lead  one  to 
anticipate  that  the  Atlantic  Water  entering  the  Mediterranean  would  be  found  along 
the  Moroccan  coast  (i.e.,  turning  to  the  right  upon  exiting  the  Strait),  in  fact 
the  inflow  is  found  near  Spain,  where  it  forms  the  northern  half  of  a  basin-wide 
anticyclonic  gyre  (Lanoix,  1974;  Cheney  and  Doblar,  1982;  Parrilla  and  Kinder, 
1984).  There  have  been  several  model  studies  of  the  gyre  (Whitehead  and  Miller, 
1979;  Nof ,  1978;  Preller  and  Hurlburt,  1982),  and  it  is  nearly  always  detectable 
in  satellite  infrared  images  (Phillipe  and  Harang,  1982). 
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OBJECTIVES 


The  primary  objective  of  the  experiment  was  to  understand  the  dynamics  and  the 
variability  of  the  anticyclonic  gyre  using  numerical  modeling,  remote  sensing,  and 
field  measurements.  Our  hypothesis  was  that  the  the  Atlantic  inflow  is  the 
primary  forcing  mechanism  for  the  gyre,  so  that  if  we  made  synoptic  measurements 
of  the  inflow  and  of  the  gyre,  and  compared  these  measurements  to  numerical  model 
simulations,  t,hen  we  could  increase  our  understanding  of  the  the  gyre  dynamics. 

The  experiment  addressed  numerous  secondary  questions  including: 

What  is  the  structure  of  the  Intermediate  and  Deep  Water  flows,  and  how  do 
these  waters  contribute  to  the  Mediterranean  outflow? 

What  is  the  biological  and  chemical  influence  on  the  gyre's  strong  ocean 
color  signal,  and  how  is  this  signal  related  to  the  physical  structure  of  the 
gyre? 


How  does  the  strong  front  along  the  northern  limb  of  the  gyre  affect  the 
atmospheric  marine  boundary  layer,  and  how  does  the  marine  boundary  layer  then 
affect  remote  sensing  data? 

How  are  the  surface  signatures  of  the  gyre,  as  measured  by  satellite  and 
aircraft  sensors,  related  to  the  subsurface  structure  of  the  gyre? 

How  does  atmospheric  forcing  influence  the  gyre,  both  directly  and  through 
changes  in  the  Atlantic  inflow? 

PRELIMINARY  RESULTS 

The  experiment  had  three  overlapping  phases,  all  ending  in  late  October  1982.  The 
first  phase  began  in  October  1981,  and  consisted  of  satellite  monitoring  (about  2 
processed  images  per  month)  of  the  thermal  surface  features.  The  second  phase 
began  in  June  1982,  when  five  current  meter  moorings  were  deployed,  a  small 
hydrographic  grid  was  occupied,  and  satellite  image  processing  increased  to  about 
once  weekly.  During  the  final  phase  in  October  1982,  four  ships  occupied 
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stations  in  the  Alboran  Sea,  the  Strait  of  Gibraltar,  and  the  Gulf  of  Cadiz  (the 
area  of  the  Atlantic  Ocean  just  west  of  the  Strait);  four  aircraft  performed 
remote  sensing  and  AXBT  flights;  a  shore-based  current-measuring  radar  (CODAR) 
operated;  additional  meteorological  and  aerosol  data  were  collected;  and  all 
available  NOAA-7  (Advanced  Very  High  Resolution  Radiometer,  AVHRR )  and  NIMBKS-7 
(Coastal  Zone  Color  Scanner,  CZCS)  images  were  obtained.  The  meeting  report  from 
an  October  1983  workshop  contains  a  detailed  account  of  the  experiment  and  some 
early  results  (Parrilla,  1984).  We  will  briefly  discuss  the  preliminary  results 
of  the  intensive  October  1982  phase. 

The  Structure  of  the  Gyre  and  The  Atlantic  Inflow.  During  October  1982  the  qyre 
was  fully  developed,  filling  most  of  the  western  Alboran  Sea.  Figure  1  shows  the 
salinity  at  100  dbar,  based  on  data  taken  over  2  weeks  from  two  ships.  The  deep 
part  of  the  inflowing  current  can  be  inferred  from  the  closely  packed  isohalines 
in  the  northwest,  and  the  gyre  as  the  low  salinity  ( <37 )  water  centered  near 
35°30'N,  4°00'W.  This  depiction  of  the  salinity  distribution  is  probably  close  to 
the  mean  for  the  period,  but  the  gyre  displayed  considerable  variability  during 
these  two  weeks.  Satellite  infrared  images  (see  cover)  showed  large  surface 
temperature  changes,  and  three  synoptic  aircraft  AXBT  surveys  showed  that  the 
thermal  center  of  the  gyre  at  100  m  moved  30  km  within  10  days. 

Near  the  Strait  of  Gibraltar,  geostrophic  estimates  of  the  Atlantic  inflow 

6  3 

relative  to  200  dbar  were  about  1.4  x  10  m  /sec  in  transport  and  maximum  speeds 
were  140  cm/ sec.  Hydrographic  stations  within  25  km  of  the  Moroccan  coast  showed 
unexpectedly  high  gradients,  and  currents  computed  from  these  data  exceeded  100 
cm/ sec. 

Current  measurements  also  showed  high  speeds  with  strong  variability  in  the  gyre 
and  the  inflow.  Five  short-period  (4-5  hour)  surface  measurements  of  sonobuoy 
drift  tracks  showed  that  the  core  of  the  inflow  was  associated  with  a  cool  sea 
surface  temperature  and  speeds  of  120  cm/sec.  South  of  the  inflow,  the  speeds  in 
the  gyre  averaged  60  cm/sec,  although  some  higher  speeds  were  measured.  Surface 
currents  measured  by  radar  (CODAR)  over  a  region  extending  60  km  south  of  Marbella 
(depending  on  interference  and  propagation  conditions)  were  well  correlated  with 
geostrophic  estimates  and  current  meter  data.  CODAR  measurements  showed  that  the 
center  of  the  inflow  current  had  a  width  of  15-30  km,  that  it  moved  more  than  30 


km  farther  offshore  in  4  days,  and  that  it  had  multiple  velocity  maxima  of  5-15 
km  width  (Janopaul  and  Frisch,  1984).  During  October,  the  current  meters  (Fig.  2) 
showed  the  subsurface  velocity  core  shifting  between  moorings  14  and  15,  with  the 
highest  speeds  at  the  shallowest  meters  (depths  of  67-124  m)  exceeding  80  cm/sec. 
The  best-resolved  velocity  section,  taken  with  the  CTD/velocity  profiler  south  of 
Estepona,  revealed  a  current  of  about  25  km  width,  100  m  depth  (20  cm/sec 
isotach),  and  surface  speeds  of  120-140  cm/sec  (Fig.  3).  This  section  also  showed 
that  in  the  highest  horizontal  shear  regions  of  the  current,  large  along-section 
velocities  made  the  non-linear  terms  in  the  momentum  equations  nearly  as  large  as 
the  Coriolis  term.  Generally,  the  geostrophic  estimates  appear  valid,  but  details 
on  scales  of  10  km  or  less  in  the  high  shear  regions  may  be  inaccurate. 

Optical ,  Biological  and  Chemical  Measurements.  Simultaneous  measurements  of 
visible  and  infrared  radiation  from  both  aircraft  (Multi-Spectral  Scanner  and 
Ocean  Color  Radiometer)  and  satellite  (C2CS  and  AVHRR)  sensors  clearly  pictured 
the  gyre  and  the  Atlantic  inflow  (Fig.  4).  The  CZCS  images  were  geometrically 
registered  to  a  Mercator  projection  and  atmospherically  corrected  for  quantitative 
chlorophyll  concentration  and  for  the  diffuse  attenuation  coefficient.  These 
results  and  the  aircraft  data  indicated  that  high  chlorophyll  concentrations  were 
correlated  with  cool  sea  surface  temperatures.  Surface  thermal  fronts  were 
coincident  with  ocean  color  fronts  on  the  large  scale,  although  differences 
appeared  at  smaller  scales.  Results  of  the  CZCS  showed  surprisingly  large  daily 
changes  in  chlorophyll  concentration  and  attentuation  coefficient  across  the 
entire  sea. 

Along  the  northern  front  formed  by  the  gyre  and  the  Atlantic  inflow,  chlorophyll 
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changed  from  1.05  to  0.45  mg/m  and  the  attenuation  coefficient  from  0.15  to  0.07 
within  24  hours.  The  cause  of  these  changes  is  still  under  scrutiny.  If  the 
cause  is  either  biological  or  physical,  then  the  rapidity,  areal  extent,  and  size 
of  the  changes  will  be  a  significant  new  result. 

In  situ  optical  and  biological  measurements  were  concentrated  along  the  mooring 
line  (Fig.  2).  These  measurements  showed  large  temporal  and  spatial  changes  in 
the  optical  and  thermal  structure.  Correlation  with  satellite  images  and  other 
data  showed  that  much  of  the  temporal  change  was  associated  with  variations  in  the 
structure  of  the  Atlantic  Inflow  and  gyre  (Arnone  and  La  Violette,  1984). 
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Biological  and  nutrient  concentrations  also  had  a  complicated  and  variable 
structure.  Phosphate  concentrations  were  higher  (0.3-0. 7  g  at  P-P04)  than 
previously  found,  and  some  high  nitrate  concentrations  also  suggested  that 
upwelling  or  mixing  may  have  been  supplying  nutrients  to  the  photic  zone.  High 
chlorophyll  concentrations  were  associated  with  the  high  nitrate.  Diatoms  were 
strongly  dominant,  with  Rhizosolenia  Stolterfothii  the  most  abundant. 

Submesoscale  Thermal  Features.  Examination  of  the  the  twice-daily  infrared  images 
(see  cover)  and  visible  CZCS  images  (Fig.  4)  revealed  tongues  of  colder  and  more 
turbid  water  that  were  pointed  toward  the  center  of  the  gyre.  They  appeared  first 
near  the  Strait,  and  then  apparently  were  advected  around  the  gyre  at  a  mean  speed 
of  40  cm/ sec  (La  Violette,  1984).  These  features  had  a  horizontal  dimension  of 
10-20  km  (compared  to  the  internal  deformation  radius  of  15-30  km),  and  a  vertical 
extent,  measured  by  AXBT  sections,  of  at  least  100  m.  Salinity  from  a  CTD  cast 
obtained  within  one  feature  was  in  the  range  of  Atlantic  Water.  During  15  days  of 
cloud-free  satellite  images,  9  features  could  be  tracked  most  of  the  way  around 
the  gyre.  Because  they  appeared  to  originate  near  the  Strait  about  twice  per  day 
(the  same  frequency  as  the  N0AA  satellite  passes),  we  conjecture  that  the  strong 
Strait  of  Gibraltar  tides  may  influence  their  formation. 

Meteorology  and  Aerosols.  Shore  stations  and  shiDboard  measurements,  including 
soundings  and  aerosols,  provided  dense  observation  of  the  lower  atmosphere.  These 
data  were  obtained  to  check  the  influence  of  the  atmosphere  on  the  gyre  (Cheney 
and  Doblar,  1982;  Bucca  and  Kinder,  1984)  and  on  remote  sensing  measurements.  The 
marine  boundary  layer  was  modified  by  cooling  and  warming  due  to  the  air-sea 
temperature  differences  over  the  water  masses  on  either  side  of  the  sea  surface 
thermal  fronts:  the  boundary  layer  over  the  colder  water  was  50  m  lower  than  over 
the  warmer  gyre  water,  had  a  stable  inversion  layer  overlaying  a  shallow  surface 
layer,  and  had  a  greater  concentration  of  smaller-sized  aerosols.  Aerosol 
populations  over  the  Alboran  Sea  (gathered  by  ship)  were  characterized  by  large 
diurnal  variations  in  concentration,  size,  and  chemical  composition.  On  the 
average,  submicron-sized  particle  concentrations  were  approximately  50%  greater 
over  the  colder  water  (probably  trapped  in  the  low  level  inversion),  while 
concentrations  of  particles  >  1  micron  were  comparable  in  both  areas.  Wave  data 
Inferred  from  aircraft  photography  also  showed  smaller  wave  heights  in  the  colder 
water.  The  correlations  of  these  atmospheric  parameters  with  oceanographic 


features  can  both  complicate  interpretation  of  remote  sensing  images  and  offer 
clues  to  useful  analysis. 

Internal  Waves.  Non-linear  internal  waves  (often  described  as  bores)  have  been 
known  in  the  Strait  of  Gibraltar  for  decades.  For  the  first  time,  however,  our 
data  show  that  these  eastward-propagating  waves  often  form  ordered  groups  of  short 
period  (about  30  minutes)  waves  in  the  Alboran  Sea,  apparently  as  internal 
solitons  (Osborne  and  Burch,  1980).  These  waves  have  sufficient  amplitude 
(exceeding  50  cm/sec  in  eastward  velocity)  to  affect  measurement  programs,  and 
perhaps  directly  influence  larger  scale  dynamics  as  well.  The  waves  were 
common,  forming  during  most  semi-diurnal  tidal  cycles  but  with  varying  strength. 

Gulf  of  Cadiz.  Hydrographic  observations  in  the  Gulf  of  Cadiz  showed  Atlantic 
Water  near  the  Spanish  coast  flowing  eastward  toward  the  Strait.  The 
Mediterranean  outflow  appeared  to  move  down  five  submarine  canyons  along  the 
Iberian  slope,  mixing  with  the  surrounding  waters.  There  were  several  gyres  or 
eddies  in  the  dynamic  topographies,  including  an  anticyclonic  gyre  (near  36°20'N 
and  7°10'W)  that  has  beer?  detected  in  CZCS  images.  Velocity  profiling  south  of 
Cape  St.  Vincent  (extreme  southwestern  Portugal)  indicated  that  the  Mediterranean 
Water  was  moving  at  speeds  of  less  than  10  cm/sec.  In  these  data  and  companion 
CTD  data,  the  temperature  and  salinity  at  the  core  of  the  Mediterranean  Water 
fluctuated  as  much  as  1°C  and  0.2  over  either  distances  of  a  few  kilometers  or 
periods  of  a  few  days  (spatial  and  temporal  fluctuations  on  these  scales  could  not 
be  separated  because  of  ship  movement). 

Numerical  Modeling.  Reduced  gravity  models  (single  active  layer  above  a  quiescent 
lower  layer)  used  by  Preller  and  Hurlburt  (1982)  demonstrated  the  importance  of 
both  the  magnitude  and  direction  of  the  inflowing  Atlantic  Water  in  the 
determination  of  the  gyre  dimensions.  Inflow  directed  north  of  east,  and  thus 
conforming  to  the  orientation  of  the  Strait  of  Gibraltar,  enlarged  the  north-south 
extent  of  the  gyre.  Increasing  the  magnitude  of  the  inflow  intensified  the  gyre 
and  shifted  it  farther  east. 

When  a  westward-flowing  lower  layer  was  Included,  the  importance  of  the  bottom 

6  3 

topography  became  apparent.  In  cases  of  large  inflows  (>1.7  x  10  m  /sec)  with  an 
Initial  northward  component,  the  submarine  ridge  system  near  Alboran  Island 


deflected  the  flow  northward.  This  deflection  eliminated  the  anticyclonic  gyre 

6  3 

west  of  the  Island.  Smaller  inflows,  such  as  the  1.4  x  10  m  /sec  estimated  for 
October,  resulted  in  a  large  anticyclonic  gyre  west  of  the  Island. 

Deep  and  Intermediate  Water.  Hydrographic  sections  showed  that  the  Deep  and  the 
Intermediate  Water  took  separate  paths  through  the  Alboran  Sea.  The  Deep  Water 
flowed  as  a  narrow  (20  km)  current  against  the  base  of  the  Moroccan  continental 
slope  and  then  appeared  at  the  eastern  end  of  the  Strait  along  the  southern  side, 
much  as  predicted  by  Bryden  and  Stommel  (1982).  The  Intermediate  Water,  however, 
appeared  to  flow  preferentially  in  the  northern  two-thirds  of  the  Alboran  Sea, 
arriving  at  the  eastern  end  of  the  Strait  along  the  northern  side.  The  Deep  Water 
flow  was  thus  anticyclonic,  like  the  shallow  Atlantic  Water  Gyre,  while  the 
Intermediate  Water  flow  was  cyclonic.  Two-layer  numerical  model  experiments 
suggested  that  this  pattern  resulted  from  the  influence  of  rotation  on  the 
Intermediate  Water,  and  of  bathymetry  on  the  Deep  Water. 
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Cover.  NOAA-7  Advanced  very  high  resolution  radiometer  (AVHRR)  infrared  images  of 
the  Alboran  Sea  collected  as  part  of  the  Donde  Va?  experiment  during  October  1982. 
The  dot  and  V  in  the  Images  show  the  advection  of  two  submesoscale  features  about 
the  Alboran  Sea  gyre  (the  line  designates  a  key  section  near  the  current  moorings, 
cf.  Fig.  2).  A  close  examination  of  the  Images  shows  that  other  unmarked  features 
were  also  being  advected.  Continuous  monitoring  of  these  cold-water  features  was 
possible  because  of  the  twice-daily  (about  every  12  hours)  spacing  of  the  NOAA-7 
overpasses.  After  registration  to  a  mercator  projection  (accuracy  +1  km)  and  an 
atmospheric  correction  to  arrive  at  absolute  temperatures,  analysis  of  the 
displacement  of  the  cold  water  features  in  successive  images  show  their  apparent 
origin  east  of  Gibraltar,  their  average  speeds  of  40  cm/sec  around  the  gyre,  and 
their  apparent  entrainment  into  the  incoming  Atlantic  Water  east  of  Gibraltar. 


1432  GMT  6  OCTOBER 


0256  GMT  7  OCTOBER 


\ 

1422  GMT  7  OCTOBER 


0244  GMT  8  OCTOBER 


1408  GMT  8  OCTOBER 


Figure  2.  Mean  velocities  from  the  current  meters.  The  solid  arrows  are  from 
instruments  in  the  Atlantic  inflow  or  gyre,  while  the  dashed  arrows  (note  seal 
change)  were  in  the  lower  part  of  the  Mediterranean  Intermediate  Water. 


Figure  3.  CTD/current  profiles  measured  south  of  Estepona  during  16-17  October. 
The  currents  were  measured,  and  not  inferred  from  the  hydrography,  (a) 

Temperature  (°C).  (b)  Salinity,  (c)  Eastward  velocity  { cm/sec),  (d)  Northward 

velocity  (cm/sec).  The  bottom  axis  is  distance  south  of  36°24'N,  and  the  top  axis 
station  location. 


Figure  4.  CZCS  image  on  7  October  1982.  Darker  shades  correspond  to  higher 
phytoplankton  pigment  concentrations  (cf.  cover). 
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Figure  5.  Pycnocline  anomaly  for  a  steady  state  two-layer  model  experiment. 

Upper  layer  inflow  was  30  cm/sec  angled  21°  north  of  east,  and  lower  layer  inflow 
was  a  uniform  0.2  cm/sec  across  the  eastern  boundary.  Positive  contours  denote  a 
thickening  of  the  upper  layer  (contour  interval  is  4  m). 
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THE  HYDROGRAPHIC  STRUCTURE  OF  THE  ALBORAN  SEA  GYRE, 
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Instituto  Espanol  de  Oceanograf ia,  Alcala  27-4°,  28014,  Madrid,  Spain 


In  1982  we  did  two  CTD  surveys  of  the  western  Alboran  Sea  as  part  of  the 
Donde  Va  experiment  (Donde  V a  Group,  1984;  Parrilla,  1984).  A  rapid  two-ship 
survey  during  5-13  October  encompassed  the  entire  western  Alboran,  and  a 
sinqle-ship  survey  during  22-30  June  covered  the  northern  half  only. 

We  use  the  hydroqraphic  distributions  and  geostroohic  calculations  to 
compare  the  qyre  found  in  1982  to  that  present  durinq  summer  1962.  Lanoix 
(1974)  used  data  gathered  by  four  ships  during  15  July-15  Auqust  1962  to 
comprehensively  describe  the  anticyclonic  gyre  which  is  formed  by  the  inflowinq 
Atlantic  Water.  The  summer  1962  and  October  1982  surveys  are  the  two  best 
realizations  of  the  qyre  to  date.  This  qyre  is  always  present,  but  its  size 
and  location  vary  (Parrilla  and  Kinder,  1984). 

We  used  two  working  definitions  for  the  gyre  boundary.  Both  the  15  dyn 
cm  (1.5  J/kq)  contour  on  the  0/200  dbar  dynamic  topoqraphy  and  the  140  dbar 
isobar  on  the  37.5  isohaline  distribution  were  closed  isolines  that  appeared 
to  separate  the  qyre  from  the  inflowing  Atlantic  Current  durinq  both  years 
( see  figures) . 

The  1962  qyre  was  larger  and  stronger  than  the  1982  gyre  (see  table). 

In  1962  the  center  was  about  14  km  farther  north,  and  the  periphery  about  20  km 
farther  north  and  30  km  farther  east.  This  resulted  in  more  than  twice  the 
surface  area  and  volume  of  Atlantic  water  (defined  as  salinity  <  37.5).  The 
1962  qyre  was  also  about  25  percent  deeper  and  had  a  stronger  dynamic  high. 
During  October  1982  transport  around  the  qyre  (i.e.,  recirculation)  was  similar 
to  the  inflow,  but  in  1962  it  was  about  50  percent  greater  than  the  inflow 
(based  on  Lanoix  Sections  8  and  9). 

The  partial  survey  in  June  1982  suggested  that  the  gyre  may  have  been 
larger  than  in  October,  althouqh  probably  smaller  than  in  1962.  The  diminution 
may  have  occurred  during  a  major  episode  of  low  eastward  flow  that  appeared  in 


ngy 
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current  meter  records  north  of  35°58'N  along  4°45'W  in  mid-September  (Kinder, 
1984).  During  this  episode  the  Atlantic  inflow  current  was  either  absent  or 
flowing  far  south  where  the  gyre  is  normally  found.  To  transform  from  a 
1962-size  gyre  to  an  October  1982-size  gyre  over  the  ten  day  low-flow  episode 
would  require  the  gyre  to  lose  about  1.1  x  10^  m^/sec  of  water  with  salinity 
<  37.5,  or  about  the  normal  eastward  transport. 
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TABLE:  Gyre  comparisons 


July-August  October  June 

1962  1982  1982  Notes 


Area  (km3) 

11,800 

4,900 

_ 

15  dyn  cm 

10,300 

4,900 

- 

37.5  at  140  dbar 

Volume  (km3) 

1,650 

725 

- 

Salinity  <  37.5 

Center:  North  • 

(35*50' 

35*40' 

_ 

dynamic  heiqht 

(35*45' 

35*39' 

- 

37.5  isohaline 

West 

(  4*08 ' 

4*10' 

- 

dynamic  heiqht 

(  4  08' 

4*10' 

- 

37.5  isohaline 

Penetration:  North 

J 36*08' 

36*03' 

36*04' 

15  dyn  cm 

(  36*08' 

35*51' 

36*05' 

37.5  at  140  dbar 

West 

(  3°.21' 

3*43' 

- 

15  dyn  cm 

\  3  18' 

3*35* 

- 

37.5  at  140  dbar 

Depth  of  37.5 

isohaline  (mean,  m) 

160 

148 

- 

(maximum,  m) 

220 

174 

>  183 

Maximum  Dynamic 

33 

25 

>  25 

Height  (dyn  cm) 

Inflow  (lO^/sec) 

1.6 

1.3 

1.2 

Gyre  (10®m3/sec) 

2.5 

1.4 

2.3 
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Summary 


Data  from  the  Coastal  Zone  Color  Scanner  (CZCS)  on  the 
Nimbus-7  satellite  have  been  processed  geometrically  and 
radiometrically  so  as  to  map  the  distribution  of  near-surface 
pigments  in  the  Alboran  Sea  during  October  1982. 


Introduction 


1  . 

During  1982  a  major  oceanographic  experiment  -  the  Alboran  Sea 
Experiment  ( "Donde  Va? ")  took  place  in  the  western  Mediterranean 
Sea  in  which  a  number  of  research  organizations  from  the  United 
States  and  from  European  countries  participated.  The  objective 
of  the  Experiment  was  the  study  of  the  jet-like  Atlantic  inflow 
into  the  Mediterranean  Sea  through  the  Strait  of  Gibraltar  and 
the  associated  anti-cyclonic  circulation  in  the  Alboran  basin. 

Since  particular  emphasis  during  the  Experiment  was  given  to  the 
verification  and  application  of  remote  sensing  methods,  relevant 
data  were  not  only  collected  from  ships  but  also  from  aircraft, 
satellites  and  ground-based  stations  in  an  attempt  to  achieve 
a  synoptic  view  of  the  complex  and  rapidly  varying  oceanic 
phenomena  in  this  region. 

Because  of  its  significance  for  investigations  concerning  both 
ocean  dynamic  processes  and  especially  marine  biology  and  ecolo¬ 
gy,  this  report  deals  with  the  remote  sensing  of  ocean  color 
from  space  and  its  evaluation  in  terms  of  the  pigment  distribut¬ 
ion  in  the  Alboran  Sea  during  October  1982.  In  Europe  these 
investigations  were  done  jointly  by  INTA  and  DFVLR. 

2 .  Ocean  color  phenomenon 
2 . 1  General 

Visible  light  may  penetrate  the  water  body  to  a  certain  depth. 
Through  processes  like  absorption,  fluorescence  and  scattering 
caused  by  natural  and  man-made  suspended  and  dissolved  substanc¬ 
es  as  well  as  by  the  water  itself,  it  is  modulated  in  intensity 
and  spectral  composition.  The  backscattered  portion  of  the  light 
as  observed  by  a  remote  sensor  above  the  water  may  thus  render 
information  on  the  type,  distribution  and  concentration  of  many 
biochemical  water  parameters.  (Because  of  the  limited^trans- 
parency  of  both  the  water  and  the  atmosphere,  only  the  spectral 
region  between  approximately  400  nm  and  700  nm  can  generally  be 
used. ) 

The  rather  complex  physical  and  biological  processes  involved 
have  been  studied  in  much  detail  during  the  past  (1).  Neverthe¬ 
less  further  research  in  this  field  is  still  needed  before  the 


interpretation  of  water  color  can  be  used  routinely  for  a 
quantitative  monitoring  of  coastal  zones.  In  particular,  the 
combined  collection  of  data  from  satellite,  aircraft  and  ship 
is  expected  to  support  the  development  and  verification  of 
algorithms  and  atmospheric  correction  methods. 

2 . 2  Optical  properties  of  water 

Radiative  transfer  theory  (2)  shows  that  the  spectral  diffuse 
reflectance  R  of  macroscopically  homogeneous,  optically  deep 
water  is  determined  at  each  wavelength  by  the  ratio  of  its 
diffuse  back  scattering  coefficient  b+  to  its  diffuse  absorption 
coefficient  a+,  as  follows: 


R 

Normally,  b+/a+<  1. 

The  spectral  reflectance  of  typical  arid  (blue)  water  results 

+  4. 

from  combined  effects  of  a  and  b  ,  since  both  vary  to  some 
extent  with  optical  wavelength  across  the  spectrum.  When  any 
type  of  living  phytoplankton  is  added  to  such  water,  the 
spectral  reflectance  rises  in  the  green  and  falls  in  the  blue 
part  of  the  spectrum.  (Compare  water-leaving  radiances  in  Fig.  1) 
The  rise  in  the  green  is  caused  by  the  increased  backscatter  of 
the  plankton  particles  increasing  b+  in  the  above  equation. 

In  the  green,  the  absorption  coefficient  a+  is  not  changed 
appreciably  by  the  addition  of  phytoplankton  because  chloro¬ 
phyll  does  not  absorb  green  light  appreciably.  Thus  the  green 
reflectance  R  increases  with  b+.  In  the  blue,  however,  chloro¬ 
phyll  absorbs  so  strongly  that  the  addition  of  phytoplankton 
causes  the  absorption  coefficient  a+  to  increase  much  more 
than  does  the  back  scattering  coefficient  b+.  Thus  R  is  decreas¬ 
ed  in  the  blue  region  of  the  spectrum.  (Compare  water-leaving 
radiance  -  dashed  line  -  in  Fig.  1) 

Sediments  in  suspension  scatter  light  strongly,  thereby  increa¬ 
sing  b+  throughout  the  spectrum,  but  they  ordinarily  have  only 
weak  effects  on  a+.  Thus  an  addition  of  sediments  tends  to  raise 
the  spectral  reflectance  of  ocean  water  at  all  wavelengths. 
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Fig .  1 


Upwelling  radiance  above  the  water  and  above  the 
atmosphere  for  clear  (solid  line)  and  for  case  I  * 
water  containing  some  chlorophyll  (dashed  line) 


Consequently  the  presence  of  suspended  sediments  can  be  disting¬ 
uished  from  chlorophyll-bearing  plankton  by  the  shape  of  their 
respective  reflectance  curves.  Decomposing  organic  materials 
produce  solutes  (yellow  substance)  which  increase  a+  primarily 
in  the  blue  and  violet  spectral  region  and  thus  lower  the  spec¬ 
tral  reflectance  at  those  wavelengths  in  a  manner  which  at  least 
theoretically  could  be  distinguished  from  the  spectral  effects 
of  living  chlorophyll. 


*  Case  I  waters  (2)  are  those  water  masses  for  which  phyto¬ 
plankton  and  their  derived  products  play  a  dominant  role  in 
determining  the  optical  properties.  They  are  not  influenced  by 
terrigenous  influx.  Case  I  waters  range  from  oligotrophic  and 
deep  blue  waters  (C  <  0.1  mg/iri3)  to  moderately  productive  green 

3 

waters  (  C^  1  mg/m  )  and  even  to  some  eutrophic,  dark  green 
waters  (  C>  10  mg/m"3)  which  can  be  encountered  in  upwelling 
regions  along  arid  coasts. 
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3  Radiative  transfer  through  the  atmosphere 


In  order  to  retrieve  the  water-leaving  radiance  (or  reflectance 
Rw)  from  the  data  collected  at  satellite  altitude,  these  have  to 
be  corrected  for  the  influences  of  the  atmosphere  in  between,  i.e. 
for  absorption  processes  and  for  scattering  by  the  molecules 
(Rayleigh  scattering)  and  the  aerosols  (Mie  scattering) . 

The  correction  scheme  is  based  on  the  following  equation  (3 )  : 


Rsat(;'>  -  +  V-'>  + 

where 


J 


(2) 


RSat(7) 

rr(^> 

V7'> 

V''> 


V7') 

ta(7') 

Tn  U) 
°3 


reflectance  measured  by  satellite 
reflectance  due  to  Rayleigh  scattering 
reflectance  due  to  aerosol  scattering 
water-leaving  reflectance 
Rayleigh  transmission  losses 
aerosol  transmission  losses 
ozone  transmission  losses 


Equation  (2)  means  that  the  reflectance  RSat(^)  at  satellite 
altitude  is  the  sum  of  two  atmospheric  terms,  R  (2)  and  R  (2)  , 
respectively  caused  by  Rayleigh  and  by  aerosol  scattering,  and 
of  the  water  reflectance  Rw(/l)  .  This  last  term  is  attenuated  by 


Rayleigh  and  aerosol  transmission  losses  TD(2) 

x\ 


and  TA ( 7 ) ;  the 


entire  term  is  in  addition  further  attenuated  by  the  factor  T^  (a  ) 

°3 

accounting  for  the  ozone  absorption. 


Further,  equation  (3)  indicates  that  the  aerosol  contribution 
RA  G* )  at  any  wavelength  2  can  be  proportionally  derived  from  a 
value  of  (Jg)  at  a  reference  wavelength  2^  , 


RA(^)  =  E(2)-Ra(2,)  .  (3) 

The  atmospheric  correction  is  the  inversion  of  equation  (2)  which 
transforms  the  satellite  reflectance  Rg  (2)  into  the  water-leav¬ 
ing  reflectance  R^(2)  by  using  equations  (2)  and  (3)  and  compu¬ 
tations  of  Rayleigh  and  aerosol  scattering  RR(*)  and  R^fa),  and 
transmissions  T  (3  )  and  Tn  (2)  . 

K  A 

This  inversion  procedure  begins  by  computing  Rc(2),  i.e.  the 
reflectance  corrected  for  ozone  and  most  Rayleigh  effects. 


Compared  to  Rgat(A)  in  which  R^x)  plays  a  major  role,  this  term 
has  the  advantage  to  be  nearly  independant  of  the  view  angles. 
Since  the  Rayleigh  and  ozone  effects  change  monotonically  across 
the  image  (in  contrast  to  aerosol  effects) ,  only  one  set  of 
conversions  from  Rgat(A)  and  Rc  CO  has  to  be  computed  and  applied 
for  all  the  pixels  of  a  given  image  sub-set.  The  size  of  this  sub 
set  has  to  be  decided  according  to  a  balance  between  correction 
accuracy  and  economy  regarding  computer  time.  For  the  images 
processed  at  DFVLR,  conversion  tables  were  computed  for  sub-sets 
of  512  x  512  pixels. 

The  introduction  of  equation  (2)  into  equation  (4)  yields  two 
further  equations: 


n>  =  Tr(2)  |r c(X)  -  Rw(^)  Ta(^)| 


Rw(A)  TA(X) 


v*> 

RC(*>  “  Td(P) 


A  Rayleigh  term  still  appears,  however,  its  effect  is  minor 

in  comparison  to  that  one  removed  in  the  previous  step.  The 
aerosol  transmittance  TA  has  a  second  order  influence.  Subsequent 
ly  R  (*)  can  be  derived  from  equation  (5)  assuming  R^^?)  is  known 
or  preferably  negligibly  small.  This  latter  assumption  is  well 
established  for  higher  wavelengths  (red  and  infrared  region) 
especially  for  case  I  waters,  and  thus  allows  the  determination 
of  R  (A) .  Once  R  (3)  has  been  determined  at  a  higher  wavelength 

A  A 

it  can  be  extrapolated  towards  the  shorter  wavelengths  (equation 
(3)  )  making  feasible  the  retrieval  of  RWU)  by  using  equation 
(6)  . 

The  factor  E()  in  equation  (3)  can  be  determined  from  the 
following  relationship: 

E  U)  =  (jr-)  0  .  (7) 

x  A(p  1 

where  n  is  the  so-called  Angstrom  exponent.  Because  n  was  not  yet 
known  at  the  time  of  the  processing,  a  value  of  n=1  was  applied. 


3 . 1  Instrument  and  orbit  specifications 

The  CZCS  on  the  Nimbus-7  satellite  is  the  first  multispectral 
image  scanner  of  high  sensitivity  dedicated  to  the  measurement 
of  ocean  color  from  space  in  five  narrow  spectral  bands. (The 
sixth  band  for  the  measurement  of  sea  surface  temperature  failed 
in  late  1981.)  The  instrument  specifications  for  the  CZCS  are 
listed  in  Table  1?  Nimbus-7  orbit  parameters  are  listed  in 
Table  2. 


FOV  . 

•  •  •  • 

•  •  • 

.  78.68  ° 

SWATH  WIDTH  . 

•  •  •  • 

•  • 

.  1566  km 

IFOV  . 

•  •  •  • 

0. 

865  m  rad 

PIXEL  SIZE  (nadir) 

•  •  • 

•  • 

. .  825  m 

SPECTRAL  BANDS: 

Centre 

Halfwidth 

Visible/near  i.r. 

443 

nm 

20  nm 

520 

nm 

20  nm 

550 

nm 

20  nm 

670 

nm 

20  nm 

750 

nm 

1 0O  nm 

Thermal  i.r. 

11.5 

F 

2  F 

Tabl.  1  CZCS  instrument  specifications 


ALTITUDE  . 

INCLINATION  . . 

.  99.2  ° 

NODE  . 

EQUAT.  CROSS  . 

. .  noon 

RE-VISIT  TIME 

ORBIT  TIME  . . . 

104  min 

ORBIT  . 

sunsynchronous 

ORBITS/DAY  . . . 

00 

• 

rn 

• 

• 

SEPARATION  . . . 

.  26  ° 

Tabl.  2  Nimbus-7  orbit  specifications 


3 . 2  Data  acquisition 

The  area  of  the  Alboran  Sea  was  covered  once  daily  by  the  CZCS 
around  noon  with  the  exception  of  4,  5,  10,  15  October  82. 

However,  not  all  of  the  imagery  collected  was  suitable  for  further 
processing.  Data  which  have  been  processed  into  chlorophyll  maps 
in  this  report  are  listed  in  Table  3. 


DATE: 

ORBIT  NR. 

APPROX.  SCAN  ANGLE 

ABOVE  ALBORAN  SEA 

1  Oct.  82 

19932 

9  °  left 

6  Oct.  82 

1  9946 

22  °  left 

7  Oct.  82 

19960 

5  °  left 

8  Oct.  82 

19974 

20  °  right 

11  Oct.  82 

20015 

35  °  left 

12  Oct.  82 

20029 

18  °  left 

1 3  Oct.  82 

20043 

6  °  right 

1 6  Oct.  82 

20084 

40  °  left 

1 9  Oct.  82 

201  26 

15°  right 

20  Oct.  82 

20140 

35  °  right 

Tabl.  3  CZCS  data  from  Oct.  82  which  have 
been  processed  at  DFVLR 


4 .  Data  processing 
4 . 1  General 

Raw  data  from  the  CZCS  are  pre-processed  at  the  Goddard  Space 
Flight  Center  according  to  NASA  standard  (Tape  Specs.  T  744041) 
into  so-called  Calibrated  Radiance  and  Temperature  Tapes  (CRTs) 
which  are  in  a  format  suitable  to  be  processed  further  on  a  digital 
image  processing  system.  The  same  standard  is  used  at  ESA's 
receiving  station  in  Lannion  and  in  Maspalomas. 

The  processing  of  CRT  tapes  was  done  jointly  by  staff  members  of 
INTA  and  of  DFVLR  at  the  Digital  Interactive  Image  Processing 
System  (DIBIAS)  of  the  Institute  for  Optoelectronics  at  the  DFVLR 
in  Oberpf af fenhofen.  For  this  purpose  a  special  software  package 
consisting  of  12  different  programmes  has  been  developed  and  in¬ 
stalled  during  the  past  years.  This  consists  of  special  read-in 
programmes,  programmes  for  additional  data  extraction,  programmes 
for  radiometric  and  geometric  manipulations  and  the  programmes 
for  the  generation  of  derived  products  like.e.g.,  chlorophyll 
maps  (4) . 

Each  CRT  contains  up  to  three  different  2  minute  scenes.  Each 
scene  consists  of  both,  image  data  (6  channels)  plus  the  addition¬ 
al  housekeeping  and  geographic  data.  These  are  essential  for  the 
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radiometric  and  geometric  processing  and  can  be  listed  or  print¬ 
ed  upon  request.  After  this,  processing  is  usually  done  on  a 
smaller,  cloud-free  sub-scene  of  the  complete  image  (approx. 

512  x  512  pixels)  and  involves  essentially  the  steps  described 
briefly  in  the  following  sections. 

4 . 2  Geometric  corrections 

The  geographic  reference  data  (anchor  points)  provided  for  each 
scan  line  are  chosen  in  such  a  way  that  they  correspond  to  the 
selected  sub-scene  (10  x  10  reference  points).  These  are  automat¬ 
ically  extracted  and  used  for  the  transformation  of  the  scene 
into  a  mercator  projection.  With  this  data  set  a  second  order 
polynomial  is  calculated  so  as  to  rectify  the  sub-scene  geometric¬ 
ally  on  a  pixel-to-pixel  base. 

Further,  a  suitable  longitude  /  latitude  grid  is  generated  with¬ 
in  the  image  processing  system  which  can  be  superimposed  onto 
the  sub-scene. 


4 . 3  Atmospheric  corrections 

In  order  to  recover  the  water- leaving  reflectance  Rw  from  the 
recorded  radiance  values  at  satellite  altitude  in  the  way 
as  described  in  section  2.3,  the  following  steps  are  used: 


o 

o 

o 

o 

o 

o 

o 

o 


Calculation  of  CZCS  sensitivity  decay  (according  to  (3)  ) 
Calibration  of  raw  data 

Conversion  into  reflectance  through  solar  constant 
Computation  of  ozone  and  Rayleigh  attenuation 
Computation  of  aerosol  reflectance  in  channel  4 
Atmospheric  correction  of  channel  4 
Atmospheric  correction  of  other  channels 
Computation  of  water-leaving  reflectances 


4 . 4  Calculation  of  pigment  distribution 

A  functional  form  which  relates  the  chlorphyll-a  concentration  C 
to  the  reflectance  at  zero  altitude  is: 


C  (mg/m3)  =  A  x  RB, 


where  R  is  the  reflectance  ratio  of  two  CZCS  channels,  while 
A  and  B  are  constants  (1). 
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For  the  processing  of  the  chlorophyll  maps  in  this  report, 
the  ratio  of  CZCS  channels  1/3  were  used  for  R,  while  for  the 
constants  A  and  B  the  values  0.47  and  -1.27  were  used  respect¬ 
ively. 

The  resulting  chlorophyll  maps  are  generated  so  as  to  cover  a 

large  dynamic  range  and  show  the  chlorophyll  distribution  in 

a  logarithmic  scale  (  white  =  high,  grey  =  medium,  black  = 

low  concentration) .  The  range  covered  in  the  imagery  corresponds 

to  chlorophyll  values  of 

0.03  mg/in  =  black 
3 

1 .00  mg/m  =  white 

The  boundary  between  land  (and  clouds)  against  water  is  masked 
with  a  white  line.  Grey  values  above  land  and  clouds  are  in 
addition  modulated  with  the  data  derived  from  channel  5,  however, 
they  are  not  significant  in  regard  to  the  information  above 
water . 

The  chlorophyll  maps  are  shown  in  Figures  2  to  11. 

5 .  Results  and  discussion 

Within  the  Alboran  Sea,  chlorophyll  appears  around  the  warm, 
anti-cyclonic  gyre,  especially  north  and  east  of  it.  This  area 
seems  to  represent  either  the  jet-like  Atlantic  inflow,  or  an 
upwelling  zone  associated  with  the  interaction  of  the  inflow 
with  the  bottom  topography  near  the  Spanish  coast  (5,  6). 

A  more  detailed  comparison  with  ship  and  aircraft  data  is  ex¬ 
pected  to  reveal  the  exact  origin  of  this  pigment  accumulation 
and  its  association  with  specific  water  masses. 

Within  the  period  during  which  CZCS  data  have  been  processed, 
rapid  changes  in  the  abundance  and  especially  in  the  geographic 
distribution  of  pigments  can  be  observed  (7) .  These  rapid  changes 
are  probably  caused  by  a  period  of  strong,  north-westerly  winds 
and  the  corresponding  displacement  and  vertical  mixing  of  water 
from  5  to  9  October  82.  Also  in  this  case,  a  more  detailed  com¬ 
parison  with  bio-chemical  (R/V  Naucrates)  and  hydrographic  data 
(R/V  Bartlett)  will  be  required  before  an  explanation  for  these 
phenomena  can  be  found. 

Generally  speaking,  the  concentration  at  which  chlorophyll 
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appears  in  the  imagery  does  not  vary  considerably  during  the 
period  of  observations.  However,  there  is  one  exception  to  this, 
namely  the  image  from  16  October  1982  (Fig.  9).  In  this  image 
no  chlorophyll  at  all  can  be  found  with  the  exception  of  the 
coastal  area  in  the  Gulf  of  Cadiz  (river  effluent) .  Since  in 
the  chlorophyll  maps  from  the  dates  before  as  well  as  from  the 
dates  after  16  October  chlorophyll  is  clearly  visible,  it  is 
very  unlikely  that  physical  or  biochemical  processes  in  the 
water  are  responsible  for  this  dramatic  change.  Instead,  it  is 
more  likely  that  the  atmospheric  conditions  were  such  that  the 
"standard"  correction  procedure  applied  to  all  data  was  not 
sufficient  any  more.  From  Table  3  it  is  evident  that  during 
the  orbit  on  1 6  October  82  the  Alboran  Sea  was  observed  at  a 
comparatively  large  scan  angle  of  almost  40  °.  Under  these 
conditions  apparently  special  steps  have  to  be  undertaken  so 
as  to  derive  corrected  water-leaving  reflectance  values  through 
the  increased  air  mass. 
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Fig.  4 


Chlorophyll-a  distribution  on 
7  Oct.  1982  (orbit  19960) 
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Fig.  5 


Chlorophyll-a  distribution  on 
8  October  1982  (orbit  19974) 
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Fig.  11  Chlorophyll-a  distribution  on 
20  October  1982  (orbit  20140) 
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THE  FLOW  OF  ATLANTIC  WATER  INTO  THE  ALBORAN  SEA  DURING  DONDE  VA 


Henry  Perkins  and  Thomas  H.  Kinder 

Naval  Ocean  Research  and  Development  Activity 
NSTL,  Mississippi  39529,  U.S.A. 

Extensive  environmental  measurements  were  made  in  the  western  Alboran  Sea 
during  5-13  October  as  part  of  the  Donde  Va?  experiment.  These  provide  a 
multi-disciplinary  and  nearly  synoptic  view  of  the  area  between  the  Strait  of 
Gibraltar  and  a  section  between  Adra  and  Cape  Tres  Forcas.  An  overview  of  the 
experiment  has  been  given  by  Parrllla  (1984)  and  by  The  Donde  Va  Group  (1984). 
Evident  in  the  dynamic  topography  is  a  large  anti-cyclonic  gyre,  similar  to  that 
reported  by  Lanoix  (1974)  but  of  lesser  Intensity  (Kinder  and  Parrllla,  1984). 

The  gyre  circulation  during  Donde  Va?  was  also  clear  in  satellite  infrared  images 
and  drifter  trajectories  (La  Violette,  1984). 

Calculation  of  geostrophic  transport  in  the  upper  200  m  (Fig.  1)  indicates  a 
satisfactory  balance  between  the  Inflow  through  the  Strait  of  Gibraltar  and  the 
outflow  to  the  east  past  Alboran  Island.  The  transport  in  the  northern  portion  of 

3 

the  gyre  along  the  section  south  of  Malaga  (2.7x106  m  /sec),  splits  into  two 
approximately  equal  parts,  one  continuing  on  towards  the  east,  the  other  returning 
westward  to  form  the  southern  half  of  the  gyre.  The  transport  at  the  eastern 
boundary  of  the  study  area  passes  in  nearly  equal  amounts  to  the  north  and  south 
of  Alboran  Island. 

Characteristics  of  the  region  in  the  northwest  where  the  gyre  and  Atlantic 
inflow  join  were  determined  by  a  series  of  stations  (Fig.  2)  made  with  an 
instrument  combining  the  features  of  a  lowered  current  meter  and  CTD  (VCTD).  Even 
at  the  westernmost  of  these  sections,  near  Estepona,  the  geostrophic  adjustment  of 
the  inflowing  water  appears  nearly  complete,  the  transport  (1.5x106  m  /sec)  of  the 
observed  current  (Fig.  3)  agreeing  well  with  that  determined  by  geostrophy. 

The  Inflowing  water  from  the  Atlantic  consists  of  North  Atlantic  Central  Water 
(NACW)  modified  by  mixing  during  Its  travel  through  the  Strait  of  Gibraltar 
(Lacomb  and  Richez,  1982).  However,  characteristics  of  the  original  NACW  can 
still  be  detected  within  the  Alboran  Basin.  Fig.  4  shows  T-S  curves  from  selected 
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VCTD  stations  along  the  Estepona  section.  Station  300  on  the  northern  edge  of  the 
current  indicates  normally  deep  water  brought  to  the  surface  through  geostrophic 
adjustment,  while  station  299  on  the  southern  edge  is  more  representative  of  the 
gyre  interior.  The  remaining  stations,  294-296,  correspond  to  the  strongest 
portion  of  the  current  (see  Fig.  3)  and  show  a  range  of  variability  expected  from 
incompletely  mixed  NACW. 

In  contrast,  T-S  curves  on  the  southern  side  of  the  gyre,  given  in  Fig.  5  for 
the  CTD  section  near  Punta  Jagerschmidt,  represent  a  well-defined  water  mass. 
Elevated  surface  temperatures  of  up  to  20.8°C  are  presumed  to  be  the  result  of 
surface  heating.  Similarly  prepared  figures  from  other  sections  show  a 
progressive  development  of  this  water  mass  as  it  sweeps  around  the  gyre. 
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Figure  1.  Distribution  of  hydrographic  stations  and  geostrophic 
transports.  Transports  are  given  in  10^  nvfysec  for  currents  be¬ 
tween  the  surface  and  the  reference  depth  of  200  m. 
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Figure  2.  Distribution  of  VCTD  stations.  Circles  indicate  ship  position 
at  the  start  of  each  station  and  a  line  indicates  its  drift  during  the 
station.  Directly  measured  current  profiles,  as  well  as  hydrographic 
parameters,  were  measured  at  these  stations. 
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Figure  3.  Section  of  eastward  current  component  for  VCTD  stations  south 
of  Estepona  during  16-17  Oct.  Station  numbers  are  given  on  the  upper  axis 
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Figure  5.  As  in  Figure  4,  but  for  CTD  stations  along  the  P.  Jager- 
schmidt  section. 
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A  TIME  SERIES  STATION  AT 
THE  EASTERN  ENTRANCE  OF  THE  STRAIT  OF  GIBRALTAR 

GREGORIO  PARRILLA 
Instituto  Espanol  de  Oceanografia 
Alcala  27,  4°  28014  Madrid,  Spain 

SUMMARY 

During  the  Donde  Va?  experiment  in  October  1982,  a  station  located  near  36°02'N, 
5°14'W  in  an  average  water  depth  of  830  m  was  occupied  for  20  hours.  Ten  casts 
were  made  with  a  continuously  recording  CTD.  The  interval  of  the  survey  included 
the  occurrence  of  two  low  tides  and  two  high  tides  at  the  Strait.  The  sampling 
rate  permitted  the  detection  of  high  frequency  oscillations  only  in  those 
occasions  when  the  down  and  up  cast  could  be  compared. 

The  isolines  showed  the  wave- like  shape,  found  by  others  (Lacombe  and  Richez, 

1982;  Cavanie,  1973),  which  is  related  to  the  tides  in  the  Strait.  A  description 
of  the  variations  of  temperature  and  salinity  with  time  is  given,  as  well  as  an 
interpretation  of  the  isolines  pattern  in  relation  with  the  flow  variations 
through  the  Strait. 

INTRODUCTION 

As  part  of  the  Donde  Va?  experiment  in  October  1982  (Donde  Va?  report,  1983; 
Parrilla,  1984),  a  time  series  of  CTD  profiles  was  made  in  the  eastern  entrance  of 
the  Strait  of  Gibraltar.  It  was  located  around  36°02’N  and  5°13'W  (Fig.  1),  close 
to  the  point  C2  of  Lacombe  and  Richez  (1982)  and  point  G  of  Cavanie  (1973).  It 
lasted  some  20  hours  (October  7,  16.00  h  to  October  8,  12.30  h)  during  a  period  of 
neap  tides.  The  average  bottom  depth  was  830  m. 

METHODS 

A  Neil  Brown  MK  III  CTD  was  used  and  calibrated  in  situ  against  water  samples. 

Due  to  the  familiar  difficulties  of  working  in  the  Strait,  the  casts  were  made 
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about  one  hour  apart,  ten  lowerings  in  total.  The  time  series  (Fig.  2a,  b)  were 
drawn  using  the  data  from  both  the  down  and  up  records  for  each  cast  (dashed 
vertical  lines).  The  time  difference  between  them  varied  from  some  60  minutes  at 
surface  to  5  minutes  at  700  m.  The  series  are  refered  to  the  tides  at  Tarifa 
(solid  vertical  lines). 

RESULTS 

The  sampling  rate  is  well  suited  to  resolve  oscillations  of  approximately  tidal 
frequency  only.  Occasionally  it  is  possible  to  notice  differences  between  the 
down  and  up  record  of  a  cast,  indicating  the  presence  of  unresolved  higher 
frequency  oscillations. 

Figures  2a  and  2b  respectively  show  variations  of  potential  temperature  and 
practical  salinity  with  time.  A  wave-like  pattern  is  seen  in  the  whole  water 
column,  although  it  is  more  uniform  in  the  upper  layer.  Its  period  conforms  to 
that  of  the  semidiurnal  tidal  component  in  the  Strait  (Lacombe  and  Richez,  1982). 

The  interface,  defined  by  salinity  37x1 0"3  and  temperature  14.5°C,  reached  its 
shallowest  value  some  4  or  5  hours  after  high  water  at  Tarifa  (HWT)  and  sank  very 
fast  for  about  1  or  2  hours  before  low  water  (LWT),  reaching  the  lowest  depth  some 
8  h  after  HWT.  The  mean  of  the  salinity  interface  oscillation  was  around  80  m. 

The  depth  of  the  interface  diminished  slowly,  taking  almost  8  h  from  trough  to 
crest,  but  it  increased  relatively  rapidly,  sinking  80  m  in  3  hours. 

In  the  upper  layers,  during  the  last  part  of  the  ascent  and  abrupt  sinking  of  the 

_3 

interface,  low  salinity  intrusions  (<36.10x10  )  and  temperature  inversions 

(Ziegenbein,  1969)  can  be  seen.  During  the  subsequent  gradual  ascent,  low 
salinities  were  not  found,  having  been  replaced  by  salinities  higher  than 

_3 

36.30x10  ,  and  the  potential  temperature  inversions  disappeared. 

In  the  lower  layer,  the  amplitude  of  the  variations  of  the  isolines  was  larger 
where  the  Levantine  Intermediate  Water  (LIW)  was  manifest.  At  casts  123  and  129, 
the  LIW  signature  was  especially  strong  between  200  and  800  meters  depth,  with 
peak  salinities  and  potential  temperatures  of  38.47x10"^  and  13.15°C  at  the  first 
of  two  extrema  and  38.45x10"^  and  13.10°C  at  the  second.  The  upper  38.43xl0~3 
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isohaline--which  together  with  its  lower  counterpart  bounds  the  LIW  core--raised 
some  200  m,  nearly  the  same  order  as  the  13°C  isotherm.  In  the  rest  of  the 

_3 

series,  the  salinity  maximum  was  always  less  than  38.44x10  and  the  depth  of  the 
isolines  was  relatively  fixed. 

However,  note  that  in  some  cases  (e.g.  casts  125  or  127)  a  variation  of  some  80  m 
has  been  observed  between  the  down  and  up  records  in  the  same  cast.  This  effect 
may  be  attributed  to  internal  waves  of  high  frequency  that  were  not  detected  in 
upper  layers  where  the  time  separation  between  both  records  was  larger. 

The  Deep  Water  (DW),  that  in  this  part  of  the  Strait  occupies  a  thin  layer  close 
to  the  bottom,  was  clearly  detected  at  cast  129  where  the  potential  temperature 
was  less  than  12.80°C  (Fig.  2a). 

DISCUSSION 

Since  there  were  no  current  measurements  simultaneous  with  the  station,  current 
data  from  other  authors  (mainly  Lacombe)  was  used  for  purposes  of  interpretation. 

According  to  Lacombe  and  Richez  (1982)  the  current  in  the  upper  layer  is  strong 
and  always  eastward,  reaching  its  maximum  near  LWT.  In  the  lower  layers,  it  is 
weaker  and  reverses  direction  with  the  tide.  The  alternating  part  of  the  current 
seems  to  be  180°  out  of  phase  between  the  two  layers.  As  a  result,  during  part  of 
the  ebb  tide  there  is  inflow  in  both  layers  reaching  a  maximum  in  the  lower.  A 
short  time  before  LWT,  the  current  in  the  lower  layer  turns  again  toward  the 
west,  while  that  in  the  upper  continues  to  increase  towards  the  east,  reaching  a 
maximum  at  about  LWT. 

During  the  flood  tide,  the  currents  in  the  upper  and  lower  layers  are  in  opposite 
directions.  After  the  outflow  reaches  its  maximum,  the  surface  inflow  decreases 
until  the  reversal  of  the  outflow  a  short  time  after  HWT,  after  which  the  water  is 
again  inflowing  in  both  layers  with  increasing  velocity. 

Taking  in  mind  this  sequence  of  currents  and  looking  at  Fig.  2,  an  intrusion  of 
North  Atlantic  Central  Water  (NACW)  (salinity  less  than  36.20x10"^)  is  seen  to 
start  at  ebb  tide  and  to  consist  of  the  Atlantic  water  which  had  accumulated 


between  this  point  and  the  western  entrance  of  the  Strait  (Ziegenbein,  1969; 

Boyce,  1975).  The  LIW  signature  was  weak  at  this  time  as  the  current  there  was 
also  inflowing.  After  the  current  reversal  in  the  ldwer  layer,  the  interface 
became  deeper,  the  most  rapid  deepening  occurring  when  the  upper  inflow  reached 
its  maximum  and  the  two  layers  system  was  no  longer  in  equilibrium  (Cavanie, 

1973).  As  the  outflow  increased,  the  LIW  was  more  evident. 

During  flood  tide,  the  inflow  in  the  upper  layer  decreased  and  the  outflow  in  the 
lower  layer  continued.  At  the  same  time,  the  interface  began  to  rise  again  and 
the  upper  layer  was  occupied  by  Atlantic  water  mixed  with  some  Mediterranean  water 
(salinity  greater  than  36.40x10”^)  and  the  layer  became  isothermal  (17°C).  When 
the  ebb  tide  began  again,  the  lower  current  became  inflowing  once  more,  initiating 
a  new  cycle. 

The  DW  was  more  visible  in  cast  129  than  in  cast  123  whereas  the  reverse  is  true 
for  the  LIW  signal.  These  differences  could  be  a  manifestation  of  the  diurnal 
tidal  component. 

The  brief  appearance  of  the  LIW  core  is  puzzling.  Why  did  its  signature  not 
persist  longer  during  the  flood  tide,  while  the  outflow  increased  and  reached  a 
maximum?  The  sampling  rate  was  barely  adequate  to  resolve  these  brief  pulses,  but 
the  fact  of  their  brevity  is  clear  enough. 
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Figure  2a.  Variation  of  potential  temperature  with  time.  Vertical  solid  line, 
high  and  low  waters  in  Tarifa.  Vertical  dashed  lines,  approximated  paths  of  the 
CTD  lowerings  and  raisings.  Horizontal  solid  lines,  depth  scale  changes. 


A  Preliminary  Study  of  a  Standing  Internal  Wave  in  the  Western  Approaches 

to  the  Strait  of  Gibraltar 


Paul  E.  La  Violette  and  Robert  A.  Arncne 
Naval  Ocean  Research  and  Development  Activity 
NSTL,  Mississippi  39529-5004 


Aircraft  flights  during  the  period  6-11  October  1984  provided  data  which 
indicate  a  standing  internal  wave  was  continuously  present  west  of  the  Strait  of 
Gibraltar  at  5°50'  West.  Although  the  wave  appeared  to  essentially  remain  at  one 
location--stretching  from  Spain  to  Morrocco,  its  surface  appearance--visual ly  and 
in  the  aircraft  radar  and  infrared  scanner  imagery— showed  considerable  changes  in 
strength  that  took  place  within  several  hours.  The  surface  feature  of  the  wave 
varied  in  appearance.  However,  it  normally  was  shown  to  have  a  main  portion  that 
consisted  of  a  central  1  kilometer-wide  ridge  of  rapidly  moving  water  with  almost 
no  waves.  On  each  side  of  this  central  region,  the  waters  were  usually  agitated 
with  a  strong  chop.  On  one  occasion,  a  wave  went  over  the  20  meter  high 
afterbridge  of  a  100  meter  long  freighter  at  the  point  that  the  freighter  left  the 
centeral  comparitavely  calm  area  of  the  standing  wave  and  entered  the  region  of 
chop. 

AXBT's  data  showed  that  horizontal  changes  of  3°C  were  present  across  the 
feature  and  that  vertical  disturbances  in  temperature  extended  at  least  to  the 
depth  of  the  AXBT's  (350  meters).  Examination  of  the  regional  bathymetric  charts 
shows  the  standing  internal  wave  was  associated  with  the  main  sill  of  the  Strait  of 
Gibraltar  and  was  evidently  a  manifestation  of  the  Mediterranean  intermediate  water 
flowing  westward  at  depth  over  the  sill  ioto  the  Atlantic. 

Unlike  the  eastward-moving  tidal-induced  internal  waves  normally  found  in  the 
strait  that  displayed  eastward-oriented  bows  as  their  surface  configuration,  the 
standing  internal  surface  wave  appearance  was  that  of  a  large  westward-oriented 
bow.  On  two  occasions,  eastward-propatating  sets  of  tidal-induced  internal  waves 
were  found  originating  from  a  position  south  of  the  center  of  the  standing  wave. 
These  appearances  occurred  approximately  2  hours  before  local  high  tide.  Although 
numerous  eastward-moving  internal  waves  were  seen  in  the  straits  during  the  flights 
none  was  noted  west  of  the  standing  wave. 
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EFFECTS  OF  WIND  VERSUS  HYDRAULIC  FORCING  ON  THE  DYNAMICS  OF  THE 

WESTERN  MEDITERRANEAN  SEA 


George  W.  Heburn 

Naval  Ocean  Research  and  Development  Activity 
NSTL,  Mississippi,  39529,  USA 


SUMMARY 


A  numerical  model  is  used  to  examine  the  relative  effects  of  wind 
versus  hydraulic  forcing  on  the  dynamics  of  the  upper  layer  circula¬ 
tion  within  the  western  Mediterranean  Sea. 


RESUME 


Les  effets  relative  du  vent  et  hydraulique  sur  la  circulation  de 
la  Mer  Mediterranee  de  l'ouest  sont  examiner  avec  un  modele  numerique. 


INTRODUCTION 


There  are  many  physical  factors  which  exert  an  influence  on  the 
circulation  dynamics  of  the  western  Mediterranean;  wind  stress, 
hydraulic  forced  flow  (inflow/outflow  through  the  Straits  of  Gibraltar 
and  Sicily),  thermohaline  circulation,  etc..  The  purpose  of  this 
study,  as  an  initial  step  in  an  investigation  of  the  hydrodynamics  of 
the  western  Mediterranean,  is  to  examine  the  relative  importance  of 
and  interaction  between  the  directly  forced  responses  to  the  two 
possibly  strongest  forcing  mechanisms,  those  being  the  wind  stress  and 
hydraulic  forced  flow.  These  two  mechanisms  are  relatively  easy  to 
study  using  a  simple  reduced  gravity,  hydrodynamic  numerical  model. 

Many  features  of  the  general  surface  circulation  of  the  western 
Mediterranean  Sea  (Allain,  I960;  Bethoux,  1980;  Ovchinnikov,  1966) 
appear  to  be  either  a  hydraulic  forced  flow,  e.g.  the  North  African 
current,  or  wind  driven  currents,  e.g.  the  circulations  in  the 
Tyhrennian,  Ligurian  and  Balearic  Seas  and  in  the  Provencal  and  north¬ 
ern  Algerian  basins.  While  it  is  true  that  the  thermodynamic  effect 
of  excess  evaporation  over  precipitation  and  river  runoff  plays  an 
important  role  in  the  Mediterranean,  this  effect  is  much  more  impor¬ 
tant  in  the  eastern  Mediterranean  than  the  western  half.  Thus  the 
western  half  of  the  Mediterranean,  divided  at  the  Strait  of  Sicily,  is 
more  amenable  to  use  of  a  purely  hydrodynamic  approach. 
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MODEL 


The  numerical  model  used  in  these  experiments  is  a  one^active 
layer,  reduced  gravity,  hydrodynamic,  primitive  equations  model  on  a 
B-plane.  The  one-active  layer  reduced  gravity  model  is  the  simplest 
model  which  can  reproduce  prominent  upper  ocean  features  and  is 
designed  to  simulate  the  first  internal  baroclinic  mode.  The  term 
reduced  gravity  is  derived  from  balance  of  vertical  forces  in  which 
the  downward  force  due  to  the  acceleration  of  gravity,  g,  is  reduced 
by  the  upward  bouyancy  force,  g(p2  -  p,)/p.  It  has  an  active  upper 
layer  separated  from  an  infinitely  deep,  inert  lower  layer  by  an 
impermeable  interface.  The  interface  represents  the  main  pycnocline 
across  which  a  density  contrast  is  specified.  Thermodynamics  and 
interfacial  mixing  are  neglected.  The  governing  equations  in  the 
vertically  integrated  mass  transport  form  in  a  right-handed  coordinate 
system  are: 
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i  and  j  are  unit  vectors  in  the  x  and  y  directions  respectively  (see 
table  1  for  a  list  of  symbols). 

The  model  equations  are  solved  using  an  explicit  version  of  the 
Hurlburt  and  Thompson  (1980)  semi-implicit  model  with  two  other  impor¬ 
tant  modifications:  1)  the  ability  to  handle  realistic  coastline 
geometry  add  by  A.  Wallcraft  (personal  communication)  and  2)  the 
outflow  boundary  condition  used  at  the  Strait  of  Sicily  is  changed  to 
a  modified  Orlanski  (1976)  radiation  boundary  condition. 


The  Orlanski  boundary  condition  is  modified  as  follows  for  use  in 
this  version  of  the  model.  The  phase  velocities  are  calculated  lo¬ 
cally  to  detect  the  presence  of  high  speed  waves  such  as  Kelvin  or 
internal  gravity  waves.  If  these  are  not  present  then  a  global  mean 
outward  advective  velocity  is  determined  using  a  mass  continuity 
integral  constraint. 

The  realistic  irregular  model  geometry  was  created  using  the 
Synthetic  Bathymetric  Profiling  System  (SYNBAPS)  (Vanwyckhouse,  1973; 
1979)  data  set  of  bottom  topography.  The  bottom  topographic  height 
values  were  interpolated  from  the  10  minute  resolution  of  the  data 
set  to  the  0.1°  by  0.05°  resolution  of  the  model. 

Table  2  lists  the  parameters  used  in  each  of  the  experiments  in 
this  study.  Based  on  observations  (Ovchinnikov,  1966;  Lanoix,  197*4; 
and  Katz,  1972),  a  200  meter  mean  at  rest  layer  depth  is  chosen  to 
represent  the  upper  layer  (layer  above  the  main  thermocline)  of  the 
western  Mediterranean.  Since  the  main  purpose  of  this  study  is  to 
examine  only  the  directly  forced  responses,  the  value  of  the  eddy 
viscosity  coefficient  was  selected  to  restrict  the  development  of 
hydrodynamic  flow  instabilities. 


FORCING  MECHANISMS 

The  model  was  driven  by  specifying  the  inflow  through  the  Strait 
of  Gibraltar,  or  by  applying  monthly  mean  wind  stresses  or  a  combina¬ 
tion  of  both.  The  inflow  velocity  was  chosen  to  yield  an  inflow 
volume  transport  on  the  order  of  1.6  Sv  (Sv  -  106  m’/sec)  . 
Estimations  (based  on  observations)  of  the  volume  transport  entering 
the  Strait  of  Gibraltar  range  from  1.2  to  1.68  Sv  (Bethoux,  1979; 
Lacombe  and  Richez,  1982;  Lanoix,  197*4)  while  the  transport  exiting 
the  Strait  of  Sicily  is  around  1.21  Sv  (Bethoux,  1979;  Garzoli  and 
Maillard,  1976).  Therefore,  considering  the  variations  of  the 
transport  estimates,  for  simplicity  and  to  prevent  the  model  from 
unrealistically  gaining  mass,  the  present  experiments  require  the 
outflow  to  match  the  inflow. 

The  wind  forcing  is  derived  from  climatogical  wind  stresses 
obtained  from  twenty  years  (1950-1970)  of  ship  observations  in  the 
Mediterranean  (May,  1982).  Individual  stresses  were  estimated  from 
ship  observations  of  wind  speed  and  direction  using  a  quadratic 
aerodynamic  drag  law  with  a  drag  coefficient  dependent  on  the  wind 
speed  and  stability.  Monthly  averages  of  the  wind  stresses  were  then 
calculated  by  averaging  the  individual  wind  stress  estimates  from  each 
month  on  a  one  degree  latitude  by  one  degree  longitude  grid.  These 
monthly  averages  were  then  bilinearly  interpolated  to  the  model  grid. 
The  wind  forcing  is  cyclic  with  a  period  of  one  year. 

The  most  prominent  and  persistent  feature  in  these  wind  stress 
fields  (fig.  1)  is  the  northwesterly  wind  blowing  through  the  gap 
between  the  Pyrenees  and  the  Alps  out  over  the  central  western 
Mediterranean.  These  wind  stresses  are  strongest  in  the  winter  and 
weaken  considerably  during  the  summer  months. 


RESULTS 


Three  experiments  were  conducted  to  study  the  relative  effects 
of  and  interactions  between  wind  stress  and  hydraulic  forcing  of  the 
Mediterranean  circulation  in  the  layer  above  the  main  thermocline.  The 
first  experiment  used  only  hydraulic  forcing  and  the  second  used  only 
wind  stress.  The  third  experiment  used  both  forcing  mechanisms. 

The  results  for  experiment  one  (fig.  2)  reveal  a  steady 
source/sink  flow  entering  at  the  Strait  of  Gibraltar  flowing  through 
the  Alboran  Sea  arcing  up  along  the  Spanish  coast  then  turning 
abruptly  eastward  through  the  Algerian  basin  along  the  north  African 
coast  and  finally  exiting  through  the  Strait  of  Sicily.  The  remainder 
of  the  western  Mediterranean  is  relatively  quiescent.  An  enlarged 
view  of  the  Algerian  basin,  i.e.  the  North  African  Current  region,  is 
presented  in  figure  3  for  later  comparison  with  the  results  from  the 
wind  only  and  combined  forcing  experiments.  This  figure  shows  the 
flow  as  it  exits  the  Alboran  with  part  of  the  flow  deflected  into  the 
Balearic  Sea  while  the  majority  of  the  flow  turns  eastward  through  the 
Algerian  basin. 

Conservation  of  absolute  vorticity  appears  to  be  the  dynamical 
driving  mechanism  for  the  path  of  the  current.  If  the  amplitude  of 
constant  absolute  vorticity  (CAV)  trajectories  (see  Haltiner  and 
Martin,  1957  pp  353-356)  were  calculated  based  on  the  current  velocity 
as  it  exits  the  Alboran  Sea  (@  35.5°N  1.0°W),  the  northernmost 
penetration  would  be  to  approximately  38°N.  The  northward  penetration 
of  the  current  (figs.  2  and  3)  is  to  around  37.5°N.  The  along  stream 
velocity  divergence  in  conjunction  with  the  relatively  high  eddy 
viscosity  could  act  to  damp  the  downstream  fluctuations  as  can  be  seen 
in  the  flow  through  the  Algerian  basin. 

The  results  from  experiment  two,  shown  in  figures  4  (full  west¬ 
ern  Mediterranean  basin)  and  5  (Algerian  basin),  reveal  quasiperiodic 
solutions  in  response  to  the  annual  cycle  of  the  wind  forcing.  The 
four  panels  of  figures  k  and  5  depict  the  depth  averaged  upper  layer 
current  velocities  for  a)  spring  (day  810),  b)  summer  (day  900),  c) 
fall  (day  990)  and  d)  winter  (day  1080).  The  variations  in  the  cir¬ 
culations  patterns  can  be  closely  related  to  the  wind  stress  curl 
variations  (fig.  6).  As  can  be  seen,  during  the  winter  when  the 
Mediterranean  is  covered  primarily  by  positive  wind  stress  curl, 
strong  cyclonic  flow  dominated  the  western  Algerian  basin,  the 
Tyhrennian,  Ligurian  and  Balearic  Seas  and  the  Gulf  of  Lions,  while  in 
the  Alboran  Sea  the  flow  is  weakly  cyclonic.  Of  particular  note  is 
the  cyclonic  eddy  in  the  western  Algerian  basin  and  the  associated 
strong  current  along  the  north  African  coast  since  this  is  in  the  same 
area  as  the  current  seen  in  the  hydraulic  forced  case.  During  the 
summer  when  most  of  the  Mediterranean  is  covered  by  negative  wind 
stress  curl,  the  cyclonic  circulations  weaken  significantly  and  in 
some  instances,  i.e.  in  the  Tyhrennian  and  Balearic  Seas  and  the  Gulf 
of  Lions,  the  flow  becomes  weakly  anticyclonic.  In  the  Alboran  Sea, 
the  circulation  becomes  moderately  anticyclonic.  Again  of  particular 
note  is  the  cyclonic  eddy  in  the  western  Algerian  basin  and 
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the  current  along  the  north  African  coast.  Both  the  eddy  and  the 
current  weaken  substantially  and  the  eddy  contracts  in  size. 

Figures  7  (full  western  Mediterranean  basin)  and  8  (Algerian 
basin)  show  the  results  from  the  third  experiment.  The  solutions  are 
essentially  the  same  as  the  wind  only  experiment  in  the  Tyhrennian, 
Ligurian  and  Provencal  basins,  while  in  the  Alboran  Sea  and  the 
Sardinian  and  Sicilian  Straits,  the  solutions  are  very  similar  to  the 
hydraulic  forcing  only  experiment. 

The  strongest  interactions  between  the  responses  due  to  the  two 
forcing  mechanisms  can  be  seen  in  the  Algerian  basin.  The  steady 
source/sink  flow  from  the  hydraulic  forcing  is  modulated  by  the  peri¬ 
odic  responses  to  the  wind  stress  forcing.  Most  notable  is  the 
interaction  between  the  hydraulically  forced  current  as  it  attempts  to 
arc  northward  along  the  Spanish  coast  and  the  wind  driven  eddy  in  the 
western  Algerian  basin.  During  the  winter  when  the  wind  driven  flow 
is  at  its  strongest,  the  flow  exiting  the  Alboran  Sea  is  forced  down 
onto  the  north  African  coast  and  the  North  African  Current  is 
intensified.  Also  note  that  the  center  of  the  eddy  is  displaced 
approximately  60  km  northeast  of  its  location  in  the  wind  forced  only 
experiment.  During  the  summer  as  the  wind  driven  flow  weakens,  the 
hydraulically  forced  current  tends  to  arc  more  northward  in  an  attempt 
to  follow  a  CAV  trajectory  both  upstream  and  downstream  of  the  eddy. 

Weaker  interactions  are  observed  in  the  Alboran  and  Balearic 
Seas.  In  the  Alboran,  the  wind  driven  flow  enhances  the  anticyclonic 
circulation  of  the  Alboran  gyre  during  the  summer  and  impedes  it 
during  the  winter.  In  the  Balearic  Sea  the  wind  driven  response 
dominates  with  a  strong  cyclonic  flow,  while  during  the  summer  the 
flow  becomes  weakly  anticyclonic  under  the  influence  of  the  hydrauli¬ 
cally  forced  flow. 


CONCLUSIONS 


A  simple  one-active  layer  reduced  gravity  numerical  model  has 
been  used  to  study  the  directly  forced  hydrodynamic  responses  due  to 
hydraulic  and  wind  forcing  on  the  western  Mediterranean  Sea.  The 
results  of  this  study  suggest  that  wind  stress  is  the  dominant  driving 
mechanism  in  the  Tyhrennian,  Ligurian  and  Provencal  basins,  while 
hydraulic  forcing  is  the  dominant  factor  in  the  Alboran  Sea  and 
Sardinian/Sicilian  Straits  region  for  the  directly  forced  solutions. 
The  flow  through  the  Algerian  basin  along  the  North  African  coast  (the 
North  African  Current)  is  significantly  affected  by  both  of  the  forc¬ 
ing  mechanisms. 

Further  investigation  into  the  hydrodynamics  of  the  Mediterranean 
will  require  the  examination  of  other  flow  regimes  within  the  physical 
parameter  space  of  the  problem.  In  particular,  the  eddy  viscosity 
must  be  reduced  in  order  to  allow  for  the  possibility  of  hydrodynamic 
flow  instabilities.  Then  a  similar  series  of  experiments  as  discussed 
here  needs  to  be  conducted  and  the  results  analyzed  to  determine  where 
and  when  flow  instabilities  are  significant. 
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List  of  symbols  for  Model  Equations 
Definition 

horizontal  eddy  viscosity 
Coriolis  parameter 
acceleration  due  to  gravity 
reduced  gravity,  g(p2  -  Pi)/p 

instantaneous  local  thickness  of  the  upper  layer 

initial  thickness  of  the  upper  layer 

time 

x-directed  component  of  current  velocity 
y-directed  component  of  current  velocity 

hv 

tangent  plane  Cartesian  coordinates:  x  positive  eastward, 
y  positive  northward,  z  positive  upward 

differential  rotation,  df/dy 

time  increment  in  the  numerical  integration 

horizontal  grid  increments 

free  surface  anomaly;  height  of  the  free  surface  above  its 
initial  uniform  elevation;  Hi  -  h  -  H 

densities  of  sea  water 

x  and  y  directed  tangential  wind  stresses  at  the  top  of  the 
upper  layer 

Table  1 
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Model  Parameters 


Parameter  Definition 


Value 


A 

eddy  viscosity 

300  m2  sec 

beta 

(df/dy) 

1.8  x  10-11  i 

f 

Coriolis  parameter 

8x10  sec 

g' 

reduced  gravity  due  to 
stratification 

-2 

0.02  m  sec 

H 

undisturbed  upper  layer  depth 

200  m 

Ax  x  Ay 

horizontal  grid  resolution 

0.1°  x  0.05° 
(7.5  x  5km) 

At 

time  step 

15  mins 

V. 

m 

inflow  velocity 

40  cm  sec  1 

Table  2 
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Figure  2. 


Upper  layer  current  vectors  for  day  1080  from  the 
hydraulic  forcing  only  experiment  for  the  full  western 
Mediterranean  basin.  Vectors  are  plotted  every  third 
grid  point  in  regions  of  relatively  strong  flow  and 
every  sixth  point  in  regions  of  relatively  weak  flow. 
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Figure  3. 


Upper  layer  current  vectors  for  day  1080  from  the 
hydraulic  forcing  only  experiment  for  the  Algerian 
basin.  Vectors  are  plotted  every  second  grid  point  in 
regions  of  relatively  strong  flow  and  every  fourth 
point  in  regions  of  relatively  weak  flow. 
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Upper  layer  current  vectors  for  a)  day  810  (spring), 
b)  day  900  (summer),  c)  990  (fall),  and  d)  day  1080 
(winter)  from  the  wind  stress  forcing  only  experiment 
for  the  full  western  Mediterranean  basin.  Vectors  are 
plotted  every  third  grid  point  in  regions  of  relatively 
strong  flow  and  every  sixth  point  in  regions  of  rela¬ 
tively  weak  flow. 


SURFACE  CURRENTS  western 


X  PL0TTED  SPEED  =  0.25  (M/SEC) 


SURFACE  CURRENTS  western 


X  PL0TTED  SPEED  =  0.21  (M/SEC) 

Fig.  4b 
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PL0TTED  SPEED  =  0. 15  (M/SEC) 


Fig. 
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Figure  5. 


Upper  layer  current  vectors  for  a)  day  810  (spring), 
b)  day  900  (summer),  c)  990  (fall),  and  d)  day  1080 
from  the  wind  stress  forcing  only  experiment  for  the 
Algerian  basin.  Vectors  are  plotted  every  second  grid 
point  in  regions  of  relatively  strong  flow  and  every 
fourth  point  in  regions  of  relatively  weak  flow. 
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Monthly  mean  wind  stress  curl  for  a)  January,  b)  April, 
c)  July,  d)  October.  Contour  interval  is 

.O 

0.2  x  10^°  dynes/ cm*. 


RANGE  IS  —  0.4E  —  06  TO  0.4E-06 


RANGE  IS  —0.1  E  — 06  TO  0.2E-06 

Fig.  6c 
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Figure  7. 


Upper  layer  current  vectors  for  a)  day  810  (spring), 
b)  day  900  (summer),  c)  990  (fall),  and  d)  day  1080 
(winter)  from  the  combined  forcing  experiment  for  the 
full  western  Mediterranean  basin.  Vectors  are  plotted 
every  third  grid  point  in  regions  of  relatively  strong 
flow  and  every  sixth  point  in  regions  of  relatively 
weak  flow. 
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Figure  8.  Upper  layer  current  vectors  for  a)  day  810  (spring), 
b)  day  900  (summer),  c)  990  (fall),  and  d)  day  1080 
from  the  combined  forcing  experiment  for  the  Algerian 
basin.  Vectors  are  plotted  every  second  grid  point  in 
regions  of  relatively  strong  flow  and  every  fourth 
point  in  regions  of  relatively  weak  flow. 
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XXIXth  Congress  and  Plenary  Assembly,  CIESM 
Lucerne,  October  1 1— 19 »  1984 
Physical  Oceanography  Committee 

Computational  Methods  for  Two  Problems  in  Air-Sea  Interaction 
B.  Edvard  McDonald,  S.  A.  Piacsek,  and  Gretchen  Dawson 
Naval  Ooean  Research  and  Development  Activity 
NSTL  Station,  MS  39529,  USA 

Abstract.  This  paper  addresses  mathematical  -  computational  techniques 
currently'  under  development  for  solving  elliptic  equations  in  the 
investigation  of  wind  -  driven  ocean  dynamics  in  large  and  small  scale 
limits. 

Resume.  Cette  communication  concerne  les  techniques  mathematiques  et 
numeriques  en  vole  de  developement  pour  la  solution  des  equations 
elliptiques  renoontrees  dans  les  investigations  de  devolution 
dynamique  des  oceans  soumis  a  l'action  des  vents  dans  la  limite  des 
grandes  et  petltes  echelles  d'espace. 

The  large  scale  problem  under  consideration 
calls  for  solution  of  an  elliptic  equation  for 
the  mass  transport  stream  funotion  T  in  a 
multilayer  rigid  lid  formulation  [1]  of 
mesoscale  wind  driven  circulation  in  a  basin  of 
irregular  shape  containing  islands.  For 

constant  depth  one  must  solve  VXV  =  c*J  in  the 
basin  with  tJ  =  0  in  the  interior  of  each 
island  Rn.  (The  constant  depth  case  is  shown 
for  purposes  of  illustration.  In  our  rigid  lid 
model [2],  topography  introduces  variable 

coefficients  into  the  equation  and  requires 
special  attention  to  layer  Intersections  with 
Boundary  conditions  are  ¥  =  0  on  the  basin  boundary^  and  =  vKn  = 
constant  on  each  island  boundary  9Rn>  corresponding  to  free  slip.  The 
values  must  be  such  that  the  circulation  around  each  island  will 
agree  with  circuit  integrals  of  the  momentum  equations.  When  applied 
formally,  this  condition  leads  to  a  capacitive  matrix  model  which,  in 
its  entirety,  is  expensive  to  use.  We  offer  an  iterative  technique  in 
which  no  circuit  integral  monitoring  is  necessary.  Our  method 
generates  sheet  vortlcity  along  each  shoreline  such  that  the  proper 
circulation  and  free  slip  conditions  are  met.  The  technique  uses  a 
symmetry  which  arises  in  a  particular  solution  for  two  dimensional  flow 
past  a  circular  cylinder:  .  . 

$  =  £  e'"1* 

where  +m  is  used  for  r<r9  and  -m  for  r>r0 .  This  solution  describes  a 
vortex  sheet  on  the  surface  r  =  r#  of  the  cylinder  with  no  flow  at 


the  sloping  bottom). 


infinity.  Vitb  properly  chosen  expansion  coefficients,  this  solution 
may  be  superimposed  upon  an  arbitrary  streamfunction  whose  vorticity 
is  zero  inside  r0  to  give  the  free  slip  condition  at  rQ ,  no  flow 
interior  to  rot  and  leave  the  flow  unchanged  at  infinity.  The  symmetry 
which  can  be  exploited  is  that  the  radial  derivative  of  S'  changes  sign 
but  not  magnitude  at  ra .  This  means  that  when  Y0  +  ¥  -  const  for  r<r0  , 
the  sheet  vorticity  at  r  =  r0  is  exactly  twice  that  which  would  be 
created  by  setting  to  zero  the  interior  flow  due  to  %.  Vith  this 
brief  justification,  one  argues  for  a  simple  Iterative  prescription  for 
the  multiple  island  problem  with  an  overrelaxation  factor  of  two.  One 
starts  with  a  solution  for  arbitrary  sheet  vorticity  at  shorelines  and 
computes  velocities  everywhere  including  land  areas.  One  then  computes 
the  vorticity  created  when  velocities  are  set  to  zero  on  land.  Exactly 
twice  this  vorticity  is  added  to  the  forcing  term  u)  and  the  iteration 
is  repeated  until  convergence  is  obtained. 

The  small  scale  problem  concerns  the  solution 
of  Poisson's  equation  for  the  velocity 
potential  for  finite  amplitude  deep  water  waves 
with  boundary  values  given  on  the  surface 
(vorticity  is  assumed  to  reside  in  a  thin  X 

surface  sheet).  A  generalized  Bernoulli 
equatlon[3]  is  used  to  forecast  the  tangential 
component  of  surface  velocity  from  known  or 
prescribed  surface  stresses.  One  solves  for 
the  velocity  potential  in  the  underlying  flow. 

This  calls  for  solution  of  a  Dlrichlet  problem  with  values  given  on  the 
arbitrarily  wavy  surface.  A  new  iteration  proceedure  allows  rapid 
solution  via  fast  Fourier  transform  after  conformal  mapping  to  a  flat 
surface.  The  mapping[4]  is  ^  +i  *%) 

X  +•  i  y  st  i  +  ^  ^ 

where  \  =  0  on  the  wavy  surface  y  =  b(x) .  The  coefficients  A  ^  are 
obtained  in  the  following  iteration  which  converges  for  wave  slopes 
less  than  unity.  }  A  ^  =  CX 

x*  =  U  * 

y*  =  ln(XH)  by  Mvrfat+ti** 

LA*  =  rPT-'  (X+iy-i)t  dl-alnseU. 
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